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Abstract 

A systematic study of the effect of fragment— fragment interaction, 
quantum statistics, 7- feeding and collective flow is made in the extrac- 
tion of the nuclear temperature from the double ratio of the isotopic 
yields in the statistical model of one-step (Prompt) multifragmentation. 
Temperature is also extracted from the isotope yield ratios generated in 
the sequential binary-decay model. Comparison of the thermodynamic 
temperature with the extracted temperatures for different isotope ratios 
show some anomaly in both models which is discussed in the context of 
experimentally measured caloric curves. 
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I. INTRODUCTION 



The response of nuclei to high excitations or temperatures has been a subject of 
intense study both theoretically and experimentally for the last several years. From 
theoretical investigations of hot nuclear matter [|l]-|3|] and also of finite nuclei [[| , it has 
been suggested that the nuclear system may undergo a liquid-gas phase transition at 
high temperatures. Recent experimental measurements of the nuclear caloric curve 
by the ALADIN Group || in the Au + Au collisions at 600 AMeV tentatively support 
such a conjecture. The key element that enters in such a surmise is the extraction of 
the nuclear temperature that they observed to be nearly constant in the excitation 
energy range of ~ 3 — 10 MeV per nucleon beyond which the caloric curve rises 
almost linearly with a slope close to that of a classical gas. Experimental data from 
the EOS collaboration ||[7| are also suggestive of critical behavior in nuclei; here too 
exact determination of the nuclear temperature has the most essential role to play. 

The temperatures of hot fragmenting systems are generally measured from the 
double ratios of isotope multiplicities employing the prescription proposed by Al- 
bergo et al [|J based on the statistical model of prompt multifragmentation (PM) ||. 
In arriving at the prescription, several simplifying assumptions are made, namely, 
(i) the fragments are non-interacting, (ii) the fragments follow Maxwell-Boltzmann 
distribution, (iii) they are formed in their ground states and (iv) all their kinetic 
energies are in the thermal mode, i.e. collective flow energy is absent. The effects 
of the interaction have later been simulated through an effective excluded volume 



interaction fL0 |; to our knowledge, the effect of fragment-fragment interaction on the 
isotope ratio temperature (T r ) within the freeze-out configuration has however not 
been taken into account. Though it is expected that at high temperatures and low 
densities the quantum system would behave like a classical Maxwell-Boltzmann sys- 
tem, the importance of invoking quantum statistics in multifragmentation has been 
emphasised by several authors [|TO| [T^|] . The qualitative effect of quantum statistics is 
to increase the number of bosons with respect to fermions at low temperatures and 
high densities, the isotope ratio and hence the extracted temperature T r might have 
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some sensitivity to the choice of statistics. The assumption of the formation of the 
fragments in their ground states is an oversimplification. In general, the fragments are 
expected to be formed in various excited states which are not too short-lived. These 
excited fragments subsequently decay either by particle or 7-ray emission. These side- 
feeding effects are shown [ pHj , P^ T5| to have an important bearing on the observed 
multiplicities and hence on the deduced nuclear temperature. The hot fragmenting 
nuclear complex that is formed in nuclear collisions may be compressed depending on 
the collision geometry which subsequently decompresses to the freeze-out configura- 
tion generating significant amount of collective nuclear flow energy. The important 
role played by collective flow on the fragmentation pattern has been shown earlier 
|l6| , [lT|j . Its effect on the nuclear temperature has only been qualitatively studied by 
Shlomo et al |18| and found to be nonnegligible. In a systematic step by step ap- 
proach, we explore in this paper the effects of the four approximations listed earlier on 
the isotopic temperatures by considering different isotope double ratios and examine 
whether they can be considered as good pointers to the thermodynamic temperature 
of the fragmenting system. 

The physics of the nuclear multifragmentation is not yet fully established beyond 
question. The one-step prompt break-up (PM) looks a very plausible scenario at 
high excitations and the sequential binary decay (SBD) model [|1^,0 may provide a 
better description of the reaction mechanism at lower excitation. Both these processes 
are thermal in nature. From the inclusive mass or charge distributions or even the 
scaling of the multiplicities of the intermediate mass fragments (IMF), it is however 
difficult [jn] to discuss the relative merits of these two competing models. If the SBD 
model is the more viable model, say, for the yield of nuclear fragments in nuclear 
collisions, then the Albergo prescription of extracting nuclear temperature from the 
double isotope ratios is called into question. One notes that in the SBD model, there 
is no unique temperature but a succession of temperatures till the nuclear fragments 
are produced in their particle stable configurations. It would still be interesting to 
know what values of temperatures one extracts from double ratios in the SBD model 
and whether they can offer some added insight in the nature of nuclear disassembly. 
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The paper is organised as follows. In Sec. II, we briefly outline the PM and SBD 
models. In Sec. Ill, temperatures calculated from both models are presented and 
discussed in the context of experimental data. The conclusions are drawn in Sec. IV. 

II. THEORETICAL FRAMEWORK 

The multiplicities of fragments produced in nuclear collisions are experimentally 
measured quantities; the nuclear temperature is a derived entity. In the following, 
we outline the models for fragment production and relate the nuclear temperature to 
the fragment yield. 

A. Prompt multifragmentation 

A hot nuclear system with N neutrons and Z protons may be formed in nuclear 
collisions at a temperature T with excitation energy e* per particle. It may be initially 
compressed in a volume smaller than its normal volume. The compressed matter 
decompresses and develops a collective radial flow in addition to thermal excitation. 
We still assume that the system evolves in thermodynamic equilibrium and undergoes 
multifragmentation after reaching the 'freeze-out' volume at a temperature T different 
from T . If the time scale involved in expansion is larger compared to the equilibration 
time in the expanding complex (i.e. the expansion is quasi-static), this assumption is 
not unjustified . We further assume that at the freeze-out volume, the system reaches 
chemical equilibrium. 

The expansion of the compressed system may be simulated through a negative ex- 
ternal pressure If there was no flow, at the freeze-out, the kinetic contribution of 
the thermal pressure is generally assumed to be cancelled by interaction contributions, 
i.e. the system is at equilibrium under zero external pressure. A positive pressure 
corresponds to compression of the system; similarly a negative pressure would cause 
decompression. If Pi is the internal partial pressure exerted by the radially outflowing 
fragments of the ith species at the surface, the total external pressure P is then given 
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by P = — J2i Pi- The total thermodynamic potential of the system at the freeze-out 
volume is given by |TB| , |2"2"| 



N 3 

G^E-TS-J^fMiUi + PQ, (1) 

i=l 

where E and S are the internal energy and entropy of the system, Q — V — Vq with 
V as the freeze-out volume and Vq the normal nuclear volume of the fragmenting 
system, N s the total number of fragment species, pi the chemical potential and u>i 
the multiplicity. The occupancy of the fragments is obtained by minimising the total 
thermodynamic potential G and is given by 

ntiPi) = exp{(e l -p l )/T}±l (2) 

where (± ) sign refers to the fermionic and bosonic nature of the fragments. The 
single particle energy is fll6l , [23" 

v 2 P 

ei = f--Bi + Vi--^. (3) 

Here E>i refers to the binding energy, pi is the density of the ith fragment species 
obtained from the momentum integration of the distribution function given by eq. (2) 
and Vj corresponds to the single particle potential , evaluated in the complementary 
fragment approximation [f24| , p5|l . It is given by 

= fexpj-U^/TpiiR^R 

Jexp[-Ui(R)/T}d 3 R 1 } 

where Ui(R) is the interaction energy of the fragment with its complementary at 

a separation R and the integration is over the whole freeze-out volume with the 

exclusion of the volume of the complementary fragment. Under chemical equilibrium, 

the chemical potential of the zth fragment species is 

Pi = p n Ni + p p Zi. (5) 

The neutron and proton chemical potentials p n and p p are obtained from the con- 
servation of baryon and charge number, iVj and Zi being the number of neutrons 
and protons in the fragment . The fragment yield is obtained from the phase-space 
integration of the occupancy function and for fermions it is given by 



^ = -^=p,X[ z Jyl fa) </>i(T) . (6) 



For bosons, the corresponding multiplicity is given by 



= 9o[e~ v ' - I]' 1 + -^n\T 3 J${r]i)<t>i{n (7) 

\/7T 1 



In eqs. (|]) and (|7|), 77^ is the fugacity defined as 



Hi + B l -V i + Pi/ pi 
Vi = f , (8) 



Aj = hj \jlTrm{r is the thermal wavelength with m 8 as the mass of the ith fragment 
species and are the Fermi and Bose integrals p6[ given by 

12 Jo ex P {(x- V )}±l W 
The first term on the right hand side of eq. (|7|) gives the number of condensed 
bosons, go being their ground state spin degeneracy. The quanity (f>i(T) is the internal 
partition function of the fragments and is defined as 

HT) = Y.9se-^ )/T (10) 

s 

where g s is the spin degeneracy of the excited state s of the cluster with excitation 



energy e*(z). The flow pressure Pi is shown to be related [|T^] to the flow energy Ei 
of the ith species in the form 

^ = C(v fi ,T)Ei (11) 

Pi 

where C is dependent on the fragment species, the temperature and also on the 
flow velocity of the fragments v^. It is found to be close to 4 except for very light 
fragments. 

In the limit 77^ << 0, (which is true when the density is very low ), J^/liv) ~^ "^ e? '> 
and then from eq. @ the yield of the fermion fragments reduces to 



Ui = nX^ 3 Af 2 eyi(T) (12) 

where At = h/\ / 2TcmT is the nucleon thermal wavelength with m as the nucleon 
mass and Ai the mass number of the ith fragment species. In the same limit, eq. (0) 



for boson yield reduces also to eq. (|T|). This is also the result obtained from the 
classical Maxwell-Boltzmann distribution. 

If one chooses two sets of fragment pairs (Ax, Zx), (A[, Z[) and (A 2 , Z 2 ), (A 2 , Z' 2 ) 
such that Z[ = Zx + p, Z' 2 = Z 2 + p, N[ = Nx + n, N' 2 = N 2 + n where n and p 
are integers, then from eq. flT^ ) it follows that the measured double ratio R 2 of the 
fragment yields can be used to determine the temperature of the fragmenting system: 

u{A! 1 ,Z' 1 )/u(A 1 ,Z l ) 



R 2 



lu(A' 2i Z' 2 )/lu(A 2 ,Z 2 ) 

(A[A 2 \ 3/2 <j>{A' 1 ,Z[,T)<l>(A 2 ,Z 2 ,T) 
\AxA' 2 J <P(Ax,Zx,T)<P(A' 2 ,Z' 2 ,T) 



e (AB/T) e -(AV/T) e (AF/T) (13) 



where 



AB = B(A[, Z[) - B(Ax, Zx) + B(A 2 , Z 2 ) - B(A' 2 , Z' 2 ) 

AV = V(A[, Z[) - V(Ax, Zx) + V(A 2 , Z 2 ) - V(A' 2 , Z' 2 ) 

AF = C[E(A[, Z[) - E(A X , Zx) + E(A 2 , Z 2 ) - E(A' 2 , Z' 2 )\. (14) 

In the limit of low density, the nuclear part of single-particle potential becomes rela- 
tively unimportant; further choosing p = and n = 1, the Coulomb contribution to 
AV practically vanishes. 

Albergo et al || further assumed the fragments to be formed in their ground states 
and they did not consider any collective flow. Then with AF = and AV = the 
temperature is easily determined from 

D / A[A 2 \ 3/2 g (A[, Z[, T)g (A 2 , Z 2 , T) 
2 \AxA> 2 ) Qq(Ax, Z x , T)g (A 2 , Z' 2 , T) 1 ) 

since the ground state degeneracy go{A, Z) and binding energies are a-priori known. 

If prompt multifragmentation is the real physical mechanism for fragment pro- 
duction, eq. (15) then provides an approximate but simple way to find out the ther- 
modynamic temperature of the disassembling system. Influences from other effects 
as already mentioned are however embedded in the experimental data for isotope 
yield ratios. One can not obtain informations on the perturbations caused by these 
effects on the double-ratio thermometer simply from the experimental isotopic yields 
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without the help of further model calculations. If there were no other effects except 
from side-feeding through 7-decay, the experimental data could be exploited to de- 
lineate side-feeding effects by using eq. (13) with AV = and AF = with the 
choice of the internal partition function from eq. (10). Effects from particle decay 
|T3| or those coming from the inclusion of Coulomb force for yield ratios involving iso- 
topes differing by proton number []27| could also be approximately reconstructed from 
the experimental fragment multiplicities. Influence of nuclear interaction, quantum 
statistics or collective expansion can not however be singled out without recourse 
to models. We have therefore done calculations in the prompt multifragmentation 
model with the barest scenario (classical statistics, no interaction, no side-feeding 
and no nuclear flow) and then included the said effects step by step to generate frag- 
ment multiplicities. The multiplicities so generated under different approximations 
are used to extract double-ratio temperatures using eq. (15) to delineate the role of 
various effects on the temperatures. 



B. Sequential binary decay 

Fragmentation may also proceed via a sequence of binary fission-like events, par- 
ticulary at relatively lower excitation energies. We employ the transition-state model 
of Swiatecki |Tj| to find the decay probability of a hot nucleus with mass A, charge 



Z and excitation energy E* into two fragments of mass and charge (Ai, Zi) and 
(A — A\, Z — Zi) respectively. At the saddle point, the binary fragmentation proba- 
bility is given by 



P(A, Z,E*;A 1 ,Z 1 ) oc exp 



2Ja{E* -V B -K)- 2VOE* 



(16) 



where a is the level density parameter taken as Aj 10 MeV' 1 , K is the relative kinetic 
energy at the saddle point and Vb the barrier height dependent on the saddle point 
temperature T s which is different from the temperature Tq of the parent nucleus given 



by To = yE*/a. The barrier height is determined in the two sphere approximation 

as 
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V B {T 8 ) = V c + V n + E sep (T , T 8 ) (17) 
where E sep is the separation energy. It is evaluated as 

E sep (T ,T s ) = B(T Q ) - 5i(T s ) - B 2 (T S ). (18) 



The binding energies are taken to be temperature dependent P5[ . The saddle-point 
temperature which is also the temperature of the fragmented daughter nuclei is given 

as 



T s = ^(E*-V B -K)/a. (19) 

The evaluation of T s from eq. ([OJ) requires a knowledge of the relative kinetic energy 
K. We assume it to follow a thermal distribution P{K) oc \f~Ke~ K l Ta . The compli- 
cated interrelationship between Vb, K and T s renders evaluation of T s difficult; to 
simplify the problem, K in eq. flI5p is replaced by its average value |T S and then T s 
is evaluated in an iterative procedure with T as the starting value. This is expected 
to be a good approximation since the dispersion in kinetic energy is of the ~ T s and 
(E* — Vb) is generally much greater than T s . The so extracted value of T s is used 
only to evaluate the barrier Vb from eq. fll7l) , the decay probability and the thermal 
distribuion. In eq. ([17]), V c is the Coulomb interaction taken to be that between two 



uniformly charged spheres and Vn is the interfragment nuclear interaction |25 . 

The relative kinetic energy K of the two fragmented nuclei lying in the range < 
K < (E*—Vb) is generated in a Monte-Carlo method obeying the thermal distribution 
as mentioned. To ensure energy conservation, this kinetic energy is plugged into 
eq. (|l^) to evaluate the temperature of the daughter nuclei for further dynamical 
evolution. The fragment kinetic energy and hence their velocities are obtained from 
momentum conservation. 

The trajectories of the fragments are calculated under the influence of Coulomb 
interaction in the overall centre of mass frame. If the fragments have sufficient exci- 
tation energy, they decay in flight. The integration of the trajectories is continued till 
the asymptotic region is reached when the interaction energy is very small (~ 1 MeV) 
and the excitation energy of the fragments are below particle emission threshold. 
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III. RESULTS AND DISCUSSIONS 



In this section we present the results of the calculations for temperatures extracted 
from double ratios of different isotope yields obtained from nuclear multifragmenta- 
tion. These calculations are performed under different approximations mentioned 
in the introduction in the PM model. For this purpose we have taken l50 Sm as a 
representative case for the fragmenting system. We also obtained the double ratio 
temperatures assuming that the fragmentation proceeds via sequential binary decay. 



A. Prompt multifragmentation 

The initial temperature T of the hot system formed is different from the kinetic 
temperature T (also referred to as the thermodynamic temperature) of the fragments 
at the freeze-out. What remains constant is the total energy E of the system or 
equivalently its excitation energy E* = E + Bq where Bq is the binding energy of the 
system. The total energy of the fragmented system may be written as 

E = ~T(M - 1) - E UiBi - i $>a>i + £ Ui < e*(i) >, (20) 
1 i=i 1 

where M = is the total number of the fragments produced in the grand canonical 

model for PM and Vj the single-particle potential. The quantity < e*(i) > is the 

average excitation energy of the ith fragment species given by 

Jepi{e)e- e / T de 
< e h) >= — — — j=-r (21) 

where the integration extend upto the particle emission threshold and is the level 



density obtained from Bethe ansatz ||20|| . To compare the temperature T and To taken 
as To = \J E* ja , we plot in fig. 1 these temperatures as a function of e* = E*/A, the 
excitation energy per particle. The dashed line corresponds to the temperature T and 
the solid and dot-dash lines correspond to the thermodynamic temperatures evaluated 
at the freeze-out volumes 6Vq and 10Vo respectively, Vq being the normal volume of 
the fragmenting system. The curve for T is parabolic but it is interesting to note that 
the caloric curves corresponding to the different freeze-out volumes mentioned show 
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plateaux in the excitation energy . In the canonical model of multifragmentation 
with multiplicity-dependent freeze-out volume, Bondorf et al [Q reported first such 
a plateau reminiscent of the onset of a phase transition in nuclei. With increase in 
freeze-out volume, we find in our calculation that the temperature decreases and the 
plateau gets extended in the excitation energy. Such a dependence of caloric curve on 
the freeze-out volume was also observed in a self-consistent Thomas- Fermi calculation 



In figures 2-7, we display the isotope double ratio temperatures T r from the prompt 
break-up of 150 Sm with different choices of isotope combinations fixing the freeze- 
out volume at 6V . The combinations are: ( 4 He/ 3 He)/(d/p), ( 4 He/ 3 He)/(t/d), 
( 7 Lt/ 6 Lt)/(d/p), ( 7 Lt/ 6 Lt)/( 4 He/ 3 He), ( 10 Be/ 9 Be)/( 4 He/ 3 He) and ( 13 C/ 12 C)/( 7 Li 
/ 6 Li). They would be referred to as (He — d), (He — t), (Li — d), (Li — He), (Be — He) 
and (C — Li) respectively. In all these figures, the dotted lines correspond to the 
temperatures obtained from the multiplicities generated in the barest Albergo pre- 
scription as mentioned earlier. It is obvious that the thermodynamic temperature 
and the double-ratio temperatures are identical in this case. The dashed lines (V int ) 
refer to the temperatures calculated from eq.(15) but with the inclusion of final state 
interaction (nuclear -l-Coulomb) over the barest scenario for the fragment generation. 
In all the cases investigated , it is found that the inclusion of fragment-fragment in- 
teraction (V) shifts the temperature by nearly a constant amount at all excitations; 
the amount of shift or its sign depends on the particular isotope combination chosen. 
The shift is found to be negligible for double ratios (He — d), (Li — He) and (Be — He). 
The dot-dash lines (QS) in the figures refer to calculations done with further inclusion 
of quantum statistics. As comparison of the dashed and dot-dash curves shows, no 
appreciable quantum effects are evident except in the case of the temperature ob- 
tained from the double ratios (Li — d). In this particular case, it is further seen that 
the difference between the quantum and classical (Maxwell-Boltzmann) calculations 
widens with excitation energy or with temperature. It is normally expected that at 



low density and high temperature fEjfl, quantum effects would not be discernible, to 



be more exact, as explained earlier it depends on whether the fugacity 77 << 0. It is 
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seen that the densities of the fragment species or alternatively their fugacity r\ vary 
in a complex way with the temperature. When the temperatue is low, the density is 
extremely low and hence the value of rj is relatively large and negative; with increase 
in temperature along with density the value of rj increases initially and then again 
decreases for the complex fragments. However for nucleons rj increases monotonically 
in the energy regime that we consider. This complex variation of r\ is reflected in 
the temperatures extracted from the double ratio of yields obtained with quantum 
statistics. 

In order to take into account effects due to side-feeding, we next assume that the 
fragments are produced in particle-stable excited states so that the ground state pop- 
ulation from the 7-decaying states have to be considered. Side-feeding from particle 



decay is thus ignored. Kolomiets et al have shown that particle-decay effects 
are rather negligible, further there is uncertainty about the cut-off limit to the par- 
ticle decay width T that one should take which is intimately coupled with the time 
scale for prompt multifragmentation. Side-feeding effects are studied after generating 
the fragment yield by using eqs.(6),(7) and (8) with flow pressure P = 0. In these 
equations, <fi is the internal partition function that includes a sum extending over the 
ground and 7-decaying excited states. For the fragments considered, isotopes upto 
4 He were taken as billiard balls with no internal excitation as it has no low-lying 7- 
decaying state. Similarly for 9 'Be, only the ground state was considered. For the rest, 
the excited states considered are 3.563 MeV for 6 Li, 0.478 MeV for 7 Li, 3.37, 5.958, 
5.960, 6.018 and 6.26 MeV for 10 Be, 4.439 for 12 C and 3.089, 3.685 and 3.854 MeV 
for 13 C. For other heavier nuclei, continuum approximation is used for the single- 
particle levels and internal partition function is taken as <p — I p(e)e~ € ^ T de where the 
integration extends upto particle emission threshold. Over and above the quantum 
statistical effects, when we consider effects due to 7-feeding, it is found from figs. 4-7 
( by comparing the dot-dash and the full lines) that these effects are very sizeable. 
The (He — d) and (He — t) thermometers show no side-feeding effects (figs. 2 and 3) 
as these fragments are taken to have no excited states. A dramatic effect is seen for 
the (Be — He) thermometer ( displayed in fig. 6) where the sharply upward going full 
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line refers to the temperature T r obtained this way. Bondorf et al Q found a similar 
behaviour for the Be — He thermometer. 

In central or near-central collisions between medium-heavy or heavy nuclei at in- 
termediate or higher energies, compression and eventual decompression of the nuclear 
matter manifests itself in nuclear collective flow energy which might be a significant 
part of the total excitation. Collective flow influences the multifragmentation pattern 



to a significant extent [IB,|IT||. The double-ratio isotope thermometer may then need 
to be recalibrated a great deal due to the nuclear flow. This is manifest from the 
figures 2-7 where the full line with crosses correspond to calculated temperatures 
with inclusion of flow above the effects induced by fragment-fragment interaction, 
quantum statistics and whereever applicable, 7-feeding. The flow energy is taken to 
be 25% of the total excitation energy. Comparison of the full line with the line with 
crosses shows that at a given excitation energy, the temperature is always lower or 
for the same temperature, the excitation energy is always higher. In fig. 8, all the 
double-ratio isotope thermometers except for Be-He are displayed for comparison. 
Except for the flow effects, other effects are included here. The behaviour of the 
temperature profiles with excitation energy look nearly the same but their magni- 
tudes differ depending on the choice of the thermometers. At lower excitations, an 
uncertainty in the T r ~ 2.0 MeV involving (Li — d) and (Li — He) thermometers is 
found which increases progressively with excitation energy. The uncertainty involv- 
ing (He — t), (He — d) and (C — Li) thermometers however decreases with excitation 
energy, all three temperatures converging at the highest excitation we study. In fig. 
9, the isotope temperature corresponding to (Li — He) is shown with inclusion of dif- 
ferent magnitudes of flow. The full and dashed curves refer to cases when half (50%) 
and one fourth (25%) of the total excitation have gone to the flow energy; the dotted 
curve corresponds to no flow. As an illustration, data from the ALADIN || and EOS 
H experiments are displayed in the figure, which use the (Li — He) and (He — t) 
thermometers respectively. To have a contact with the EOS data, we also display 
the calculated temperature from the (He — t) thermometer with 50% flow energy 
(dot-dash curve). In an analysis of the same data in Ref [29fl, it was pointed out that 
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the data could be better explained invoking progressive increase of the percentage of 
flow energy with increasing total excitation; comparison of the present calculations 
with the experimental data validates this observation. 

B. Sequential binary decay 

Hot nuclear systems may release energy through binary fission-like decay, the 
decay chain continues till there is no further energy for binary division. At the end 
of such decay process, fragments of different species are produced in ground states 
and in 7-decaying excited states, the multiplicity depending on the initial system and 
excitation energy. It has been noted earlier [[3(J that the frequency distribution of the 
fragments follows almost a power-law distribution and that it is not too different from 
the one obtained from prompt multifragmentation at the same excitation energy. Our 
calculations done at different excitation energies also show that the inclusive mass or 
charge distrbutions obtained from both PM and SBD models are roughly the same. 
The isotopic distributions are however seen to have significant differences. In the SBD 
model, the hot nucleus prepared initially at an excitation energy or temperature goes 
through a succession of decays, the temperature of the produced fragments (assuming 
equilibration before each decay) therefore also decreases as time proceeds. In fig. 10, 
we display the average temperature T av of the produced fragments as a function of 
time when the initial system 150 Sm has been prepared at three different excitation 
energies, namely, e* = 13.5, 10.0 and 6.0 MeV. The temperature of the fragments 
is calculated from T av = (10 < e* >) 1//2 where < e* > is the ensemble averaged 
excitation energy per particle of the fragments at any particular instant of time. It 
is found that the higher the initial excitation energy of the system, the faster is the 
cooling rate which is expected. An experimentalist does not know a-priori whether 
multifragmentation is a one-step process (PM) or is an outcome of a sequence of 
binary decays. If one takes the fragmentation yields from the SBD model as the 
'experimental data' , it would be interesting to see the results for the double ratio 
temperatures calculated with the Albergo prescription as given by eq. (15). The 
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double ratio temperatures so calculated for the combinations (He — d), (He — t) 
(Li — d) and (Li — He) are displayed in fig. 11. One finds that except for (Li — d), 
the temperatures are very weakly dependent on the initial excitation energy and are 
very low (~ 3 MeV) even at the highest excitation energy we study. Such apparent 
temperatures were obtained by Ma et al in their Albergo-type analysis of the 
experimental data in 36 Ar+ 58 Ni collisions at 95AMeV. For the (Li—d) thermometer, 
the temperature however rises steadily with initial excitation. Thus the functional 
dependence of the temperature T r with excitation energy obtained from the SBD and 
PM models are very different; the thermometers in the SBD model also register too 
low a temperature compared to the PM model. 

The kinetic energy distribution of the fragments at the end of the decay process 
would reflect the overall kinetic temperature of the system. In the SBD model, 
since the system proceeds through a sequence of temperatures, the kinetic energy 
distribution reflects an apparent temperature. In fig. 12, this apparent temperature 
Tkin is shown as a function of initial excitation energy from the slope of the final energy 
distributions of p, d, 3 He and A He produced from 150 Sm. The temperatures extracted 
from the four distributions are not very different. Closer inspection however shows 
that except for the one for A He, the 'Caloric curves' show broad plateaux mimicing 
a liquid-gas phase transition. This arises possibly from the changing temperature 
scenario and a complicated energy dependence of the fragment partial widths for 
decay in the SBD model. 

IV. CONCLUSIONS 

We have calculated the apparent temperatures from several combinations of the 
double-ratio of isotope yields in two different physical scenarios in perspective; the 
one-step prompt multifragmentation and the sequential binary decay. In the PM 
model, the inclusion of final state interaction gives rise to nearly a constant shift 
in the temperature T r calculated as a function of excitation energy from the one 
obtained from the Albergo prescription, the shift being different for different isotope 
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combinations. The effect of quantum statistics on the apparent temperatures is found 
to be nominal; the effect of 7-feeding is very substantial and is found to be rather 
dramatic for the (Be — He) thermometer. The presence of collective flow reduces 
the apparent temperature T r for a given total excitation energy. Moreover, a soft 
plateau, generally seen in the caloric curves obtained for the double-ratio temperatures 
becomes extended with inclusion of flow energy. The import of our calculations is 
that better contact with the experimental data can be achieved if one assumes that 
the excitation energy has a collective flow component in it. 

One can not rule out the sequential binary decay as a possible reaction mechanism 
for the fragment yields, particularly at not too high excitation. This prompted us 
to study the caloric curves where the apparent temperatures T r are calculated from 
the fragment yields in the SBD model, both from the double-ratios and slopes of 
the energy distributions of the fragments. The double ratio temperatures generally 
show extended plateaux but no subsequent rise at higher excitations; on the other 
hand the caloric curves calculated from the slopes of the energy distributions display 
broad shoulders with subsequent rise at higher excitations mimicing a first-order phase 
transition. Since caloric curves obtained in both the PM model and the SBD model 
show apparent signatures of a phase transition, conclusion regarding phase transition 
in nuclear collision requires utmost caution and search for additional signatures is 
called for. 
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Figure Captions 

Fig. 1 The temperature of the fragmenting system 150 Sm as a function of e*, the exci- 
tation energy per nucleon. The dashed line (T ) corresponds to the temperature 
in the Fermi-gas approximation, the full and dot-dash lines refer to temperatures 
at the freeze-out volumes taken to be 6V and 10V respectively. 

Fig. 2 The temperature T r from the double-ratio ( A He/ 3 He)/(d/p) for the system 
150 Sm as a function of e* in the prompt multifragmentation model. The dot- 
ted line refers to T r obtained from Albergo prescription , the dashed, dot-dash 
and full line with crosses correspond to the temperatures with subsequent pro- 
gressive inclusion of final state interaction, quantum statistics and flow energy, 
taken to be one-fourth of the total excitation. 

Fig. 3 Same as figure 2 for the ( 4 He/ 3 He)/(t/d) thermometer. 

Fig. 4 The temperature T r from the yield ratio ( 7 Li/ 6 Li)/(d/p) in the PM model. 
The dotted line refers to the Albergo prescription, the dashed, dot-dash, full 
line and line with crosses refer to T r with subsequent step by step inclusion of 
final state interaction, quantum statistics, 7-feeding and flow energy, taken to 
be one-fourth of the total excitation. 

Fig. 5 Same as figure 4 for the ( 7 Li/ 6 Li)/( 4 He/ 3 He) thermometer. 

Fig. 6 Same as figure 4 for the ( 10 Be/ 9 Be)/( 4 He/ 3 He) thermometer. 

Fig. 7 Same as figure 4 for the ( 13 C / 12 C) / ( 7 Li/ 6 Li) thermometer. The dashed and 
dot-dash curves can not be distinguished from each other. 

Fig. 8 The double-ratio temperatures T r obtained after inclusion of final state inter- 
action, quantum statistics and 7-feeding in the PM model. The fragmenting 
system is 150 Sm. The solid, dashed, dotted, dot-dash and the full line with 
open squares refer to (He — d), (He — t), (Li — d), (Li — He) and (C — Li) 
thermometers respectively. 
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Fig. 9 The temperatures calculated with all effects as discussed in the text from the 
(Li — He) and (He—t) thermometers in the PM model. The fragmenting system 
is 150 Sm. The dotted, dashed and full lines refer to calculations for (Li — He) 
temperatures with 0%, 25% and 50% of the total excitation as flow energy. The 
dot-dash line refers to the (He — t) temperatures with flow energy as 50% of 
the total excitation. The crosses refer to ALADIN data and the open squares 
refer to the data from EOS collaboration. 

Fig. 10 Evolution of the average temperature T av as a function of time for the 150 Sm nu- 
cleus prepared at excitations of 13.5, 10.0 and 6.0 MeV per nucleon respectively 
in the SBD model. 

Fig. 11 The double-ratio temperatures from the thermometers (He — d), (He — t), 
(Li — d) and (Li — He) when the fragments have been produced in the SBD 
model from 150 Sm. 

Fig. 12 The kinetic temperatures obtained from the energy distributions of the frag- 
ments p, d, 3 He and 4 He produced in the disassembly of 150 Sm in the SBD 
model. 
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